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Abstract 

A thin film superconductor shaped into a spherical shell at whose center lies 
the end of long thin solenoid in which there is an integer flux iV<I>o has been 
previously extensively studied numerically as a model of a two-dimensional 
superconductor. The emergent flux from the solenoid produces a radial B- 
field at the superconducting shell and N vortices in the superconducting film. 
We study here the effects of including a second solenoid (carrying a flux /) 
which is inserted inside the first solenoid but passing right across the sphere. 
This Aharonov-Bohm (AB) flux does not have to be quantized to make the 
order parameter single valued. The Ginzburg-Landau (GL) free energy is 
minimized at fixed TV as a function of / and it is found that the minimum is 
usually achieved when the AB flux / is half a flux quantum, but depending 
on N the minimum may be at / = or values which are not obvious rational 
fractions. 

PACS numbers: 74.20.De, 74.60.-w 
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A spherical shell with a magnetic monopole at its center provides a useful geometry for 
numerical studies of a thin film (two-dimensional) superconductor in a perpendicular mag- 
netic field!. Furthermore, an investigation of the ground state of the vortices in this system 
revealed an interesting geometric effect. While the ground state of vortices penetrating an 
infinite type II superconductor plane is the well-known Abrikosov vortex lattice!, where 
vortices form a triangular array, on a spherical surface a perfect triangular lattice cannot 
form without the presence of at least twelve disclination defects and when N, the number 
of vortices on the surface is large other defects appear!. 

In this paper, we descibe a situation where again geometry and topology play key roles. 
We consider as in Ref. [I] the ground state of vortices in a spherical superconductor in the 
presence of a radial magnetic field generated by a monopole at the center of the sphere, 
but the vector potential describing the magnetic field contains also an additional Aharonov- 
Bohm (AB) flux. It should be possible to realize this system experimentally by inserting 
a solenoid into the center of a spherical superconducting shell as the field which emerges 
from the end of the solenoid approximates to the field from a magnetic monopole. One may 
visualize this system as a spherical superconductor with a couple of very thin solenoids one 
of which ends at the center of the sphere and the other lies along the z-axis. (See Fig. [I].) 

The quantum mechanical system which consists of a magnetic monopole and an AB flux 
is known to have many unusual properties as discussed in Ref. [|in detail. For instance, there 
exist solutions to the Schrodinger equation for which the Dirac quantization of monopole 
charges! does not hold even if the wavefunction is required to be single-valued. In our 
system, the monopole charges should be quantized, since it is determined by the number of 
vortices penetrating the superconductor as discussed below. One of the main results of this 
paper is that, in the ground state, the system organizes itself such that a nonvanishing AB 
flux is induced with its strength given by a fraction of the fundamental flux quantum with 
the actual value related to N, the number of vortices, in a way which is obscure to us. This 
'quantization' of the AB flux differs in its origin from other types of flux quantizations which 
are usually determined from a topological consideration such as the requirement that the 
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order parameter remain single valued as one passes round a circuit. In the present system, 
the flux is quantized dynamically in the sense that it is determined by minimizing the free 
energy for a superconductor. 

We consider a spherical type II superconductor of radius R and width d <C R in the 
presence of a radial magnetic field H = H(r)r. This magnetic field is produced by a 
magnetic monopole at the center of the sphere. The system can be described by the following 
Ginzburg-Landau free energy for a (complex) superconducting order parameter ty(9,(j)): 



where dVL is the solid angle element, a,f3 and m are phenomenological parameters, and 
D = — iV — (e*/hc)A with the vector potential A and the charge of a Cooper pair e* = 2e. 
Physical properties of the system including the effect of fluctuations are described by the 
partition function given by 



We study the ground state of the system by employing the mean-field theory in which one 
approximates the partition function in @ by the saddle point values of F, and neglects all 
the fluctuations. Furthermore we shall work in the extreme Type II superconducting limit. 
For the vector potential (or for the magnetic induction B = V x A), this amounts to taking 
the magnetic induction such that its value is equal to the magnetic field from the monopole 
at the superconducting surface. In terms of the vector potential, one could take the following 
standard form for a magnetic monopole with charge g, which has a Dirac singularity along 
the negative z-axis: A r = Aq = 0, = (g/r) tan(#/2). Note that B = g/R 2 . However, a 
more general form for A will be considered in the functional integral in (0) by adding an 
Aharonov-Bohm potential which describes the effects of an infinitely thin solenoid carrying 
flux / along the z-axis in addition to the magnetic field produced by the monopole. On the 
surface of the superconductor, the vector potential in this case can be written as 




+ ^|^|4 + ^|vx A-H| 2 1, 
2 1 1 8tt 1 J ' 



(1) 




(2) 
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A r = A e = 0, = BRtan - + J , (3) 

2 ZttK sin U 

Note that except on the z-axis this vector potential gives the same B as in the case without 
the AB flux. The presence of such AB potential, however, makes a profound change in the 
spectrum of the operator D 2 in a spherical geometry^. 

As for the order parameter, we employ the lowest Landau level (LLL) approximation, 
which one usually uses at the mean-field theory level to describe the vortex lattice!. We 
expand \I/ in terms of the eigenstates of D 2 and discard all the higher eigenf unctions except 
the lowest one. For a system consisting of N vortices, the flux quantization condition, which 
results from the single-valuedness of the order parameter around a vortex, gives that the 
total flux through the surface of the sphere, B(4ttR 2 ) is equal to N times the fundamental 
flux quantum $o — hc/e*. This condition is equivalent to the Dirac's quantization condition 
for the monopole charge g: 2e*g = Nhc. The normalized LLL wavefunctions are given byi 

W,<M = ^e^sinl^^^cos^-™-' 1 ^), (4) 

where an integer m labels the degeneracy of the LLL, a = —f/<&o, and the normalization 
constant h m = [AnR 2 B(\m + a\ + 1, \N — m — a\ + l)] -1 ^ 2 with the beta function B. The 
corresponding eigenvalues of D 2 are [1(1 + 1) — N 2 /A]/R 2 , where 

/ = -[\ m + a \ + \N - m - a\]. (5) 

Since zeros of the order parameter are identified with vortices, we might expect that iV+1 
eigenfunctions are needed to describe the positions of N vortices and an overall constant in 
the order parameter. We note then that, because of the peculiar form of the eigenvalues 
given in ([5]) , one has to include at least one eigenfunction which has a larger eigenvalue than 
the rest of the eigenfunctions. Let us first focus on the case where < a < 1. The remaining 
cases will be discussed later. We may take eigenfunctions with m = 0, 1, 2 ■ • • iV (which is 
not a unique choice as discussed below) to form a LLL order parameter as follows: 

N 

f,^p^^ m (y ;ffi ), (6) 

m=0 



where P = {ksThc/ dfiBe*) 1 ^ and v m are complex numbers. Inserting this into ([!]) we 
obtain 

-^[{^m}] r v-^ I |2 . / \| |2l 

7 rp = <*t[ 2^ \ v m\ +e{a)\v N \\ 

m=0 

1 W 

S k 1+ k 2 ,k 3 +k 4 W 12 34V* kl vl 2 V k3 V k4 , (7) 



2iV fcl =o 



where 



dP 2 . he*B^ 

CUT = -, — 7^(Oi H ) 

k B T y 2mc' 

is the dimensionless temperature which changes sign at the mean field transition line, H c2 (T), 

2a, 1 + a, dP 2 he*B , s 

e (a) = 1 + — 1 + — - — -- , 8 

a.T iV kb± 2mc 



and 



_ g(^Eifci + 2a + l^Eil^-fci-a| + l) (q] 
Wl2U IlMh + a +l,\N-ki-a\ + 1)]V2 1 j 



Finally, if we rescale t> m = |aT| 1//2 w m , we have 



-^[{^m}] 2 r I 12 /Ml 

. rp = a T[ - 2^ \ u m\ -e{a)\u N \ 
m=0 



1 N 



5 k 1+ k 2 ,k 3 +k4 W 1234U* kl U* k2 U k3 U k4 \, (10) 



2 ^ ki =o 

where we are mainly concerned with the region below the mean field transition line (i.e. 
ax < 0) as the minus signs on the right hand side indicate. Within mean field theory the 
vortex lattice configuration is determined by minimizing the free energy (|10"D for given N. 
We denote the minimum free energy per vortex (F /NhsTa^) for given a and N by En (a). 
Once the coefficients {u m } that minimizes the free energy are obtained, the positions of the 
vortices can be determined as elaborated later. 

We first discuss an interesting symmetry relation in this system. We note that in (||) one 
could also use ip m with m = — 1, 0, 1 • • • N — 1 for which one eigenf unction (with m = — 1) 
has a larger eigenvalue than the rest. If we start from 



N-l 

i>{6,(f>) = P v m ij> m (e,tia), (11) 

m=— 1 

one can easily show that 

N 

i;(9, 0) = S N -WP v m i> m {n ~ 0, -0; 1 - a), (12) 

m=0 

where v m = v^-m-i- Inserting this into (|]) we find that the free energy expression is just 
the same as (|7|) except that a is replaced by 1 — a. Therefore the minimum free energy 
En (a) as a function of a has a reflection symmetry about a = 1/2 for < a < 1: 

E N (l-a)=E N (a). (13) 

Eq. (|l^) also suggests that the vortex configuration at 1 — a can be obtained by performing 
the transformations; 9 — > 7r — 8,<fi — > —0 on the vortex lattice obtained at a. 

There is also a periodicity relation with respect to a in the following sense. For 1 < a < 2 
we have two choices: either m — — 1, 0, 1 - • -N — 1 or m = —2, — 1, ■ ■ ■ iV — 2. If we take 
the first and start from (|TT|) , we can show that 

v 

0) = e-*P £ « m Vm(0, 0; « - 1), (14) 

m=0 

where t> m = t> m _i in this case. Proceeding as in the previous case, we conclude that 

E N (a) = E N (a - 1), (15) 

and the vortex lattice has the same configuration for both cases. One can easily see that 
([T5p holds for all a (including integral values of a). Together with ( |T3"D it characterizes the 
ground state energy spectrum as a function of the flux strength. Keeping these relations in 
mind, we will restrict our discussion below only to the region < a < |. 

We study the ground state of the vortex system described by ([]]) and (|||) by numerical 
minimization of (fT0|) . We use a straightforward quasi-Newtonian algorithm. The free energy 
(p!0|) is invariant under the transformations, u m — > exp(i7)w m , and u m exp(imj)u m , for 
arbitrary real 7, where the latter amounts to the rotation of the vortices by an angle 7 
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around the z-axis. These symmetries enable one to set, for example, Imuo = Imui = to 
reduce the number of independent variables from 2N + 2 (N + 1 complex variables) to 2N. 

The first thing we notice in the numerical minimization is that the minimum free energy 
occurs when vn = 0, i.e., the mode that gives the larger eigenvalue than the rest is not 
present after all. This means that, for < a < 1, we can write the order parameter as 

^(MHPsin^cos 1 -^) 

N-1 a a 

x £ v m h m e^sm m {-) C os N ~ l - m {-) 

m=0 1 1 

= Csin a (^)cos 1 - a (^) (16) 

N-1 a a 

x II {a k M^ /2 -b k cosC-)e-^ 2 ), 
k=i 1 1 

where a k = cos(8 k /2) exp(-z</> fc /2) and b k = sin(0 fc /2) exp(i0 fc /2) with (8,(f>) = (8 k ,4> k ) 

being the positions of vortices. The overall complex constant C and N — 1 pairs of real 

numbers (8 k , <p k ) account for the original iV complex coefficients v m . Note that the presence 

of the AB flux drives the vortices into a peculiar configuration described by flT6|), where in 

addition to iV — 1 vortices at (8 k , <p k ) there are fractional vortices at the north (6 = 0) and 

south (8 = tt) poles in the sense that the zeros are of order a and 1 — a respectively. In 

the limit where a — > we obtain a configuration where one full vortex at the south pole. 

As already noticed in the previous studies0 for a = 0, one is free to fix the position of one 

vortex, since the free energy depends only on the relative positions of vortices. Therefore the 

continuation from the nonzero a to a = uniquely picks the positions of vortices among the 

degenerate configurations. When one continues the order parameter to a = 1, one reaches 

the configuration where one full vortex is on the north pole. This configuration is related 

to the one for a = via (|I2"D, and is among the many degenerate ground states of vortices 

without an AB flux. 

We obtain the ground state energy -Eat (a) as a function of a for various values of N. As 
discussed above, it is sufficient to consider the interval, < a < |. The main results of our 
investigation is that for most of the values of iV we studied the ground state energy E^(a) 



has a minimum at nonzero values of a. This means that the vortices organize themselves 
such that a flux of strength / = — a$o is induced along the z-axis. Typical examples of 
energy spectra are shown in Fig. |2|. For iV = 1,2, • ■ -36, we have found two exceptional 
cases at iV = 14 and N = 22 for which the energy spectrum also shows a local minimum 
or maximum as shown in Fig. |3|. While we have not made a systematic investigation of 
the scale of the variation of En (a) with AB flux / as N increases, we suspect that it is 
largely N independent. If so then the overall free energy change scales as N, which for 
large N could dwarf the magnetic energy stored in the AB solenoid and justify our use of 
the word "induced". In other words the AB flux in this limit could be regarded as being 
spontaneously generated. 

For some values of N, the ground state energy E^(a) has a minimum at a = i.e at 
vanishing AB flux. These cases include N = 2,3, 6, 12, 32, 42 • • •. We note that except for 
the first three, these numbers correspond to the so-called magic numbers found in previous 
studies on the vortex system on a sphere^ without the AB flux. It was found there that at 
these numbers the ground state energy has a lower value than for nearby values of N, and 
the vortex configuration in this ground state possesses a five-fold symmetry. 

The nonzero value of the induced AB flux in most cases is equal to one half of the 
fundamental flux quantum (a = |). In particular, for all odd values of N up to N = 35, the 
free energy is minimized for a — |. But, for some values of N, the value of a is given by 
interesting fractions other than ~. These cases up to iV = 36 are tabulated in Table. |. 

The actual ground state configuration of vortices when the AB flux is present is also of 
interest. We obtain the configuration of vortices that minimizes the free energy for given 
value of a by numerically finding the roots (9k, (pk) from the coefficients {u m }. We study 
how it changes when the flux strength a is increased. As remarked above, for a ^ there 
are fractional vortices at the north and south poles with strength a and 1 — a respectively. 
We find in general that, since the strength of the fractional vortex at the north pole becomes 
increased as the value of a increases, the nearby vortices are pushed downwards. We expect 
that for most cases one reaches a kind of balanced configuration between vortices on the 

8 



northern and southern hemispheres as one approaches a = ~, which contributes to lowering 
the energy. An example of this behavior can be seen in Fig. |j. A word of caution should be 
noted at this stage, however, that the movement of vortices as the value of a is changed is 
not entirely a well-defined concept, since the vortex configuration is only determined up to 
certain symmetry relations: A rotation around the z-axis by any amount or a reflection in 
the x-z plane (which corresponds to u m — > w^J does not change the energy. In Fig. |, for 
instance, in order to compare the vortex configuration for two different values of a, we have 
rotated the vortices until the three on the equator in both cases coincide. 

Even if one fixes the positions vortices in an appropriate way as in Fig. |], the detailed 
movement of vortices as the AB flux strength is increased is not a simple downward transla- 
tion but more complicated. This is especially true for the case where the minimum energy 
is achieved for nonzero fractions other than | as listed in Table [p. In this case, while the 
fractional vortices on the poles give the nearby vortices the tendency to move downwards 
(or upwards), the other vortices also organize themselves in a complicated way to achieve 
the ground state. An example of this behavior can be seen in Fig. [5]. We believe there must 
be a topological reason (yet to be understood) which relates these nonzero fractional values 
of the AB flux strength to the particular values of N. 

In summary, we have presented an interesting example of an interplay between the ge- 
ometry of a superconductor and the vortices on it resulting in the quantization of an AB flux 
with its strength given by a fraction of the fundamental flux quantum. We stress again that 
this quantization of an AB flux has a dynamical origin following from the minimization of 
the free energy for a vortex system. Since our analysis is based on the mean field theory, our 
results will be valid only in the very low temperature regime. As the temperature is raised, 
we expect the quantization effect will be weakened by thermal fluctuations and the AB flux 
will fluctuate around its minimum value. It would also be interesting to study possible 
relations of the present results to other physical systems exhibiting fractional vorticesil. 
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TABLES 

TABLE I. Fractional values (other than 1/2) of the Aharobov-Bohm flux quanta up to A = 36. 



A a = -f/if> 

8 0.1305 ± 0.0005 

14 0.2775 ± 0.0005 

20 0.4365 ± 0.0005 

22 0.1895 ± 0.0005 

24 0.1915 ± 0.0005 

30 0.3250 ± 0.0005 

34 0.3545 ± 0.0005 

36 0.3615 ± 0.0005 
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FIGURES 




FIG. 1. A spherical superconductor penetrated by two solenoids: (b) provides an approximately 
radial magnetic field at the spherical surface and (a) carries an AB flux along the z-axis. 
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FIG. 2. The ground state energy per vortex as a function of the AB flux strength a = —f/Q>o 
for various values of N; (a) N = 12, (b) N = 24, (c) N = 28, (d) N = 36. The dotted lines are 
guides to the eye 
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FIG. 3. The ground state energy per vortex as a function of the AB flux strength for (a) 
= 14, and (b) N = 22. 
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FIG. 4. The configuration of vortices for N = 28 and for a = 0.05 (filled circles) and for a = 0.5 
(open circles). The case where a = 0.5 has a lower energy. Figures (a) and (b) describe the vortices 
on the northern and southern hemispheres respectively. The poles correspond to the centers of the 
circles, on which fractional vortices are represented by the smaller size of the symbols. 
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FIG. 5. The configuration of vortices for TV = 36 and for a = 0.15 (filled circles), a = 0.25 (open 
triangles), a = 0.35 (open squares), and a = 0.45 (open circles). The case where a = 0.35 has the 
lowest energy. To compare, we have rotated the vortices until the nearest one to the north pole 
lies on the x-axis. (a) and (b) represent the vortices on the northern and southern hemispheres as 
in Fig. g. 
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